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Abstract 

We present a lattice QCD calculation of the p meson decay width via the P-wave scattering phase 
shift for the 1 = 1 two-pion system. Our calculation uses full QCD gauge configurations for A'^^ = 2 
flavors generated using a renormalization group improved gauge action and an improved Wilson 
fermion action on a 12'^ x 24 lattice at m^^/mp = 0.41 and the lattice spacing 1/a = 0.92 GeV. 
The phase shift calculated with the use of the finite size formula for the two-pion system in the 
moving frame shows a behavior consistent with the existence of a resonance at a mass close to the 
vector meson mass obtained in spectroscopy. The decay width estimated from the phase shift is 
consistent with the experiment, when the quark mass is scaled to the realistic value. 
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I. INTRODUCTION 



A study of the p meson decay is a significant step for understanding the dynamical as- 
pect of hadron reactions with lattice QCD. We notice three studies carried out to date 
toward this direction. The earlier two [l|, [2I employed the quenched approximation ignor- 
ing the decay into two ghost pions which appear in the quenched theory. The third one 
used the full QCD |3j and estimated the decay width from the p ^ tttt transition ampli- 
tude (pItttt) extracted from the time behavior of the correlation functions (0|7r(t)7r()f:)p(0)|0) 
and (0|7r(t)7r(t)7r(0)7r(0)|0), assuming that the hadron interaction is small. These studies, 
however, were all carried out at unphysical kinematics m^/mp > 1/2. 

In the present work we attempt to carry out a more realistic calculation. We estimate 
the p meson decay width by calculating the P-wave scattering phase shift for the 1 = 1 
two-pion system. The calculations are carried out with Nf = 2 full QCD configurations 
previously generated for a study of the light hadron spectrum with a renormalization group 
improved gauge action and a clover fermion action aX j3 = 1.8, k, = 0.14705 on a 12^ x 24 
lattice j^. The lattice parameters were determined from the spectrum analysis which gave 
mT^/nip = 0.41, the lattice extent L = 2.53 fm and the lattice space inverse 1/a = 0.92 GeV. 
The finite size formula presented by Rummukairrerr and Gottli;b i is employed to esfmate 
the phase shift. This calculation is made at two energies that allow to study the existence 
of the resonance. 

This paper is organized as follows. In Sec. [TTlwe give the method of the calculations and 
the simulation parameters. We present our results in Sec. IIIII Our conclusions are given 
in Sec. IIVI Preliminary reports of the present work were presented in j^. The calculation 
was carried out on VPP5000/80 at the Academic Computing and Communications Center 
of University of Tsukuba. 

II. METHODS 

A. Rummukainen-Gottlieb formula 

Let us consider the p meson decay into the two pions in the P-wave. When the p meson 
is at rest, the energy of the two pions, neglecting the final state interaction, is 



where the momentum takes p = (27r/L)n (n G 7^ 0) on the lattice. In a typical full QCD 
simulation on L ~ 2 — 3 fm, this energy is significantly larger than the resonance mass rrip. 
On our full QCD configurations, for example, the lowest energy estimated from m^r and nip 
calculated in the previous study is P = 1.47 x mp, appreciably away from the resonance 
and it is not suitable to study the p meson decay. 

In order to realize the kinematics such that the energy of the two pions is close to rrip, 
we consider a system having a non-zero total momentum, i.e., the moving frame [3], with 
the total momentum p = pe^ = [271 / L)e^ in a box satisfying the periodic boundary 




(1) 
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condition. We set the system to the representation of the rotation group on the lattice 
(the tetragonal rotation group I^ih), which represents the J = 1 spin state ignoring effects 
from higher spin states with J > 3. We consider the iso-spin representation of (/, J^) = (1, 0), 
the neutral p meson. 

For non-interacting hadrons the dominant low energy states in the moving frame are the 
two free pions with the momenta p and 0, and the p meson with the momentum p and a 
polarization vector parallel to the momentum, P3(p). The energies of these states are 



Wi = Y/m^~+~p^ + ttIt, for the two free pions , 

W2 = ^/iTi'j^ + p2 for the p meson . (2) 

Other states having higher energies are neglected. On our full QCD configurations the 
invariant mass of the two free pions takes y/s = 0.97 x nip, which is closer to rup than that 
given by ([T]) for the system having the zero total momentum. 

The hadron interaction shifts the energy from to Wn {n = 1,2), and the energies 
Wn are related to the two-pion scattering phase shift 6 in the infinite volume through the 
Rummukainen-Gottlieb formula 0], which is an extension of the Liischer formula j3] to the 
moving frame. The formula for the representation and the total momentum p = pe^ 
reads 

-— - = Z(l;A;L/(2vr)), (3) 
tano 

where k is the momentum defined from the invariant mass -^i as ^/s = — = 

2 a/^^ + . The function Z is an analytic continuation of 

^5:i±(Mi. (4, 

which is defined for Re(a;) > 5/2, where 7 = W/y/s is the Lorentz boost factor and the 
summation for r runs over the set 

• I n = ni , r2 = ^2 , rg = + /7 , n G | . (5) 



Z{l;q) can be evaluated by the method described in Ref. [81]. 



B. Extraction of energies 

In order to calculate the two energies Wn (n = 1,2) we construct a matrix of the time 
correlation function, 

= ( (0| {^'^)\t) {ir'n){ts) |0) (0| (7r7r)t(t) p,{ts) |0) \ 
^ ^ V (0| p\it) i7rn){ts) |0) (0| pl{t) p,{ts) |0) J ' 

where Psit) is an interpolating operator for the neutral p meson with the momentum p = 
{2n/L)e3 and the polarization vector parallel to p, and (7r7r)(t) is an interpolating operator 
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for the two free pions, 



(7r7r)(t) = i=(7r-(p,t)7r+(0,t)-7r+(p,t)7r-(0,t)) . (7) 

These operators belong to the and the (/, J^) = (1, 0). 
To extract Wn {n = 1, 2) we construct a matrix, 

M{t,tR) = G{t)G-\tR) , (8) 

with some reference time tn j^. The two eigenvalues Xn{t,tR) {n = 1,2) of the matrix 
M{t, tji) behave as 

A„(t,tK)=e-^" (*-*«) , (9) 

for a large t, if the two lowest states dominate the correlation function. The two energies 
Wn can be extracted by a single exponential fit to A„(t,t/j). 

In order to construct the meson state with a non-zero momentum we introduce U{1) 
noise ^j(x) in three-dimensional space, which satisfies 

Nr 

J2C*{^)Cj{y)=S%^-y) for Nn ^ oo , (10) 

i=i 

where is the number of the noise representation taken to be 10 in the present work. We 
calculate the quark propagator, 

QABMq,ts,Q = J2iD-')AB{^,t;y,ts) ■ [e*''%(y)] , (11) 

y 

where A and B refer to color and spin indices. The square bracket in ffTTl) is taken as the 
source term in solving the propagator. The two point function of the meson with the spin 
content F and the momentum p can be constructed from Q by 

E E^"'''" (75Q^(x,t|o,t5,0)75rtQ(x,t|p,t5,0)r) , (12) 

j=l X 

where the bracket means the trace with respect to the color and the spin indices; hereafter 
we take this convention for the brackets. 

The contractions of the quark field for the components of G{t) are shown in Fig. [H The 
vertices refer to the pion or the p meson with the momentum specified in the diagrams. The 
time runs upward in the diagrams. The function G^jr^jr7r(^) for the first diagram in Fig. [T]is 
calculated by introducing in addition another U{1) noise ^^^(x) having the property identical 
to that of ^j(x) as in (fTOj) . 

Nr 



i=l x,y 



(13) 
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The function G^7r^7r7r(^) for the second diagram is obtained by exchanging the momenta of 
the sink in (fT^ . 

To obtain the G{t) for the other diagrams we calculate the quark propagator of a different 
type by the source method, 

WABi^,t\k,t,\q,ts,ij) = J2J2^D-')Aci^,t-z,ti) ■ Q{z,tM,ts,^M , (14) 



z C 



where A, B and C refer to color and spin indices, and the term in the square bracket is 
taken as the source term in solving the propagator. Using W we calculate the functions 
(^^^^^^^^(t) for the third to sixth diagrams in Fig. [T]by 



^[3rd] _ 


Nr 

j = l X 


-ipx 


(W\^, t\0, ts\ - p, ts, QW{^, t\0, t\0, ts, 


(7 [4th] _ 


Nr 


-ipx 


(w{^, t\0, ts\p, ts, QW\^, t\0, t\0, ts, q) 


(7 [5th] ^ 

TTTT — ^TTTT 


Nr 

EE^- 

J = l X 


-ipx 


(w{^, t\p, ts\o, ts, Qw^{^, t\o, t\o, ts, 


^[6th] _ 

TTTT — ^TTTT 


EE'- 

j = l X 


-jpx 


(W\^, t\ - p, ts\0, ts, QW{^, t\0, t\0, ts, 


tions Ctttt- 


for the two diagrams in Fig. [T]can be calculated by 


r<[ist] _ 


Nr 

EE^- 

j = l X 


-ip-x ^ 


^Q(x, t\0, ts, ^j)W\^, t\ - p, t^lO, ts, ei)7573 


n[2nd] _ 

TTTT— >p 


Nr 

EE^- 

j = l X 


-ip-x ^ 


^iy(x, t\p, ts\0, ts, QQ\^, t\0, ts, , 



J/ / 5 



(15) 



(16) 



and similarly for Gp^nirit), 



Nr 

Gi^^L = -EE^"*"" (w\^,t\0,t\0,ts,QQMp,ts,^,Hi3] 

j = l X 

Nr 

= -^^e-P-'^ (^Q^M-p,ts,^j)WMO,t\0,ts,Qi,is) . (17) 

j=l X 

The quark propagators are solved with the Dirichlet boundary condition imposed in the 
time direction and the source operator is set at = 4 which is sufficiently large to avoid 
effects from the temporal boundary. We calculate the Q-type propagators ( ITTl) for four sets 
of q and the U{1) noise: 

(q, noise) = { (0,0, iO,v), (P,0, (-P,0 } • (18) 
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The W-tjpe propagators (fT^ are calculated for 22 sets of k, ti and q: 

(k,ti|q) = { (p,ts|0), (-P,t5|0), (0,t5|p), (0,t5|-p), (0,ti = 4-21|0) } , (19) 

using the same U{1) noise representation ^ in common. All time correlation functions can be 
calculated with combinations of these propagators. We carry out additional measurements 
to reduce statistical errors using the source operator located at ts + T/2 with the Dirichlet 
boundary condition at T/2, and average over the two measurements. Thus we calculate 
(4 + 22) X 10 X 2 = 520 quark propagators for each configuration. 

C. Simulation parameters 

Calculations employ Nj = 2 full QCD configurations previously generated for a study of 
the light hadron spectrum using a renormalization group improved gauge action and a clover 
fermion action at /5 = 1.8, k, = 0.14705 with the mean-field improvement taking Csw = 1-60 
on a 12^ X 24 lattice |5]. The periodic boundary conditions are imposed for both spatial and 
temporal directions in configuration generations and the Dirichlet boundary condition for the 
temporal direction in calculations of quark propagators. The lattice parameters determined 
from the spectrum analysis are mT^/mp = 0.41, L = 2.53 fm and 1/a = 0.92 GeV. The total 
number of configurations analyzed every 5 trajectories is 800. We estimate the statistical 
errors by the jackknife method with bins of 100 trajectories. 

III. RESULTS 

A. Time correlation function 

In Fig.[2]we show the real part of the diagonal components (tttt — > nn and p ^ p) and the 
imaginary part of the off-diagonal components (tttt — >■ p and p tttt) of the time correlation 
function G{t) in ([6]). G(t) is a Hermitian matrix, since the sink and source operators are 
identical for a sufficiently large Nr or equivalently for a large number of configurations. The 
off-diagonal components are pure imaginary by P and CP symmetry. We find that these 
hold true within statistics. The p ^ Tin component agrees with mi —>■ p a.s seen in Fig. [2] 
within the error, but the statistical errors of the former is large for a large t. Hence, in the 
following analysis we substitute p tttt by tttt — p to reduce errors. 

We calculate the two eigenvalues Xnit,tR) {n = 1,2) for the matrix M{t,tji) in (IHl) with 
the reference time = 9. In Fig. [3] we plot the normalized eigenvalues Rnit, tji) {n = 1, 2) 
defined by 

R^itM) - >^nit,tn)-^-^—^^-^ , (20) 
where G.,t{t; p) is the time correlation function for the pion with the momentum p, 

a(t;p) = (0|7rt(p,t)7r(p,ts)|0) . (21) 
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The slope of the curves in Fig. [3] represents the energy difference with respect to the energy 
of the two free pions, i.e., = Wn — W^. We observe that the energy difference for 

RiitjtR) is negative and that for R2{t,tji) positive, meaning that the phase shift is positive 
and negative, respectively, consistent with the presence of a resonance in between. 

We extract the energy difference AWn for both states by a single exponential fit to 
Rn(t, tn) for the time range t = 10 — 16. The energy of the two free pions is calculated 
from the mass m,r and the energy E obtained by a single exponential fit to GTr{t;0) and 
G^{t; p) in ([2ID, as = + E. The energy Wn is reconstructed by = AWn + W^. 
The results are tabulated in the upper part of Table [H 

B. Effect of finite lattice spacing 

To see the size of errors arising from the discretization of the energy and the momentum 
on the lattice, we study the accuracy of the dispersion relation for the single particle. In 
the upper part of Table [III we show the mass m and the energy E of the pion and the p 
meson with the momentum p = (27r/L)e3, calculated from the time correlation functions. 
The p meson with zero momentum Pj{0) (j = 1, 2, 3) and that having the momentum with 
the perpendicular polarization Pj{p) (j = 1,2) cannot decay energetically, so that the mass 
and the energy are obtained from the time correlation functions in a usual way. 

In the continuum we have a relation 

E = + ^2 ^ (22) 

for the single particle. We expect that the relation is modified on the lattice to 

cosh(E) = 2 ■ sin^(p/2) + cosh(m) . (23) 

For the measured m and E we calculate the momenta Pes from the two dispersion relations 
fl22]) and fl25]) . We expect that these p^s agree with the pre-set momentum p = 27r/L up to 
the discretization error. The lower part of Table HTl shows Pcs and the ratio pl^/p'^ whose 
departure from unity is a measure for the violation of the dispersion relation due to the 
discretization error. We find that the relation on the lattice (l23l) satisfies well, whereas that 
in the continuum fl22|) does not hold so well for the pion. For the p meson the departure 
from unity is not clearly detected due to a large statistical error. 

We should also be concerned with the discretization error for the two-pion system. We 
may think of two sources of the errors in the Rummukainen-Gottlieb formula ([3]). One of 
them arises from the Lorentz transformation from the moving frame to the center of mass 
frame using Lorentz symmetry in the continuum. In the transformation we use the relations, 

= ^W^ - p2 , 
e = s/4 - ml , (24) 

for the invariant mass ^/s, the energy in the moving frame W and the momentum k. These 
relations suffer from the discretization error, and hence the definitions of y/s and k contain 
similar errors. 
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After the Lorentz transformation we obtain the Helmholtz equation for the two-pion wave 
function 0(x) in the center of mass frame, 

(V^ + k'^^(j){x) = for|f|>i?, (25) 

where R is the two-pion interaction range. In Ref. j3] Rummukainen and Gottheb derived the 
formula by solving fl25|) ignoring the effect of the finite lattice spacing. Thus the discretization 
error also appears here. 

The violation of the dispersion relation in the continuum (122|) has also been shown in 
Ref. 01, where the phase shift for the S'-wave state was calculated for a statistical model. 
Motivated from the validity of the dispersion relation on the lattice (|23|) . Rummukainen 
and Gottlieb obtained the invariant mass ^/s and the momentum k from the energy in the 
moving frame W as 

cos h(v^) = cosh(iy) - 2 ■ sin^(p/2) , 
2 ■ sin2(A;/2) = cosh(v^/2) - cosh(m^) , (26) 

and the phase shift was obtained by substituting k into their formula. The discretization 
error arising in solving (!25|) . however, was not worked out. 

In the present work we calculate the invariant mass ^/s and the momentum k from the 
energy-momentum relations both in the continuum (IMl) and on the lattice (l26ll . and estimate 
the phase shift by putting k into ([3]). We take the difference arising from the two choices as 
the discretization error. It is expected that this error vanishes in the continuum limit. 

C. Scattering phase shift and decay width 

The invariant mass ^/s, the momentum k and the phase shift S obtained by the procedure 
in the previous subsection are presented in the lower part of Table [B Appreciable differences 
that depend on the energy-momentum relations used are visible in ^/s and k, but the 
difference for 6 is comparable with statistical errors. These are also shown in the lower 
panel of Fig. |U where the phase shift sin^ 6, which is proportional to the scattering cross 
section of the two-pion system, is plotted. In Table [T] we see that the sign of 5 at < mp 
{nip = 0.858 ± 0.012) is positive (attractive interaction) and that at ^/s > rup is negative 
(repulsive interaction) as expected. This confirms the existence of a resonance at a mass 
around nip. 

It may in principle be a straightforward task to estimate the p meson decay width by 
fitting the phase shift data with the Breit-Wigner formula. The quark mass we worked with, 
however, is much heavier than the realistic value, so that a long extrapolation is needed. 
Since we expect kinematic factors in the decay width depend largely on the quark mass while 
we made a simulation for only one set of the quark mass, we avoid this direct measurement 
of the decay width and take a different approach. We parametrize the resonant behavior of 
the P-wave phase shift in terms of the effective p ^ tttt coupling constant ^fp^^r as 

g'^TTTT k^ 

tan 5 = -^—j^TT-ps r , (27) 
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as in the continuum theory, where Mr is the resonance mass and gpT^jr is defined by 

LeS = gpnn eabcih - k2)f, pl{p)'K\ki)'K\k2) . (28) 



abc 



We may expect that such a couphng constant does not vary too rapidly as the quark mass 
changes. 

In (l27j) the invariant mass \/s and the momentum k that appear in the effective theory 
in the continuum satisfy the energy-momentum relations in the continuum fl2^ . A further 
discretization error may then appear in the definition of and k upon the application of 
(!27|) to the phase shift calculated on the lattice. We find, however, that this does not lead 
numerically to a serious error. In Tabled we give the momentum k^ = s/4 — calculated 
from ^/s, and show that the difference between k and fco is small. Thus the error caused by 
the different definitions of the momentum is negligible. We adopt ko in the application of 

(EZD. 

The results of the coupling ^fp^^ and the resonance mass Mr given by fl27j) are 

Qp^^ = 6.25 ±0.67 
Mr = 0.851 ±0.024 
Mr/tUp = 0.992 ±0.033 (29) 



using the energy- momentum relations in the continuum (]24l) . and 

gp^^ = 5.82 ±0.55 
Mr = 0.906 ±0.028 
Mr/tUp = 1.056 ±0.038 (30) 

using those on the lattice fl26|) . where rrip is the p meson mass extracted from the time 
correlation function in the previous section. In the lower panel of Fig. H] we draw the curve 
for sin^ 6 given by (1271) with gp-^^, and Mr given in ( l29l) and ( l30l) . The position at sin^ 6 = 1 
that corresponds to the resonance mass Mr is also plotted in the upper panel of Fig. IHfor 
the two cases and compared with nip. We find that Mr is consistent with rrip. 

Assuming that the dependence of gpj^T^ on quark mass is small, we estimate the p meson 
decay width at the physical quark mass as 

pph ^ 9%n (fcPy ^ 2 ^ 4_^28 MeV, (31) 

where = 770 MeV is the actual p meson mass and (k^^)^ = (mP'^)^/4 — (mP^)^ {m^ = 
140 MeV). This yields 

FP^ = 162 ± 35 MeV (32) 



for ([29]), and 

pph ^ ^ 27 MeV (33) 

for (1301) . These estimates are consistent with experiment, F = 150 MeV. The difference 
that arises from the two energy-momentum relations used is comparable with the statistical 
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error. This is an encouraging result, although we assumed that the coupling constant does 
not depend on the quark mass to make a long extrapolation and we did not trace the 
propagation of discretization errors that appear in various steps so accurately. 

IV. CONCLUSIONS 

We have shown that a calculation of the P-wave scattering phase shift for the 1=1 
two-pion system and estimation of the decay width therefrom are feasible on the lattice 
with present computing resources. The phase shift data shows the existence of a resonance 
at a mass close to the vector meson mass obtained in the spectroscopy. This resonance can 
certainly be identified with the p meson. We extracted the p meson decay width from the 
phase shift data and showed that it is consistent with the experiment. 

We note, however, at the same time several important issues that should be cleared in 
the future work. One of the important issues is to reduce the discretization error in the 
kinematic relations among the invariant mass ^/s, the energy in the moving frame W and 
the momentum k on the lattice, needed to evaluate the phase shift. An obvious way to solve 
this problem is to use a lattice closer to the continuum limit. 

We have used the effective p — > tttt coupling constant gpT,-,, to extrapolate from the point 
mT^/nip = 0.41, where our simulation is made, to the physical point rriT^/mp = 0.18, assuming 
that (jfpTTTr does not depend on the quark mass. More direct evaluation of the decay width is 
clearly desirable. The decay width may be estimated directly from the energy dependence 
of the phase shift data by fitting the Breit-Wigner resonance formula, if the simulations 
are made close to the physical quark mass and we have data for several energy near the 
resonance mass. We must leave these issues to studies in the future. 
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FIGURES 
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FIG. 1: Quark contractions of vrvr vrvr, vrvr p and p vrvr components of the time correlation 
function G{t). Vertices refer to the pion or the p meson with a momentum specified in the diagram. 
The time runs upward in the diagrams. 
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FIG. 2: Real part of the diagonal components (vrvr tttt and p ^ p) and the imaginary part of 
the off-diagonal components (vrvr p and p — > tttt) of the time correlation function G{t). 
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FIG. 3: Normalized eigenvalues Ri{t,tR) and R2{t,tR). 
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( upper panel ). Cont refer to the results obtained with the energy-momentum relations in the 
continuum (j24p and Lat to those with the relations on the lattice (j26p . The two lines are given by 
([27|) with parameters gpj^.,^ and Mr given in ([29j) and ([30]) . The abscissa is in lattice units. 
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TABLES 





n = 1 




n = 2 




" n 
Wn 


0.9805(51) 
-0.0441(60) 
0.9364(63) 




0.105(22) 
1.085(22) 




V~s 

tan 5 
sin^ 5 


Cont 

0.7764(76) 
0.0235(26) 

0.07906(93) 
0.00621(14) 


Lat 

0.8000(77) 
0.0337(28) 
0.0328(27) 
0.0736(66) 
0.00539(96) 


Cont 

0.951(25) 
0.099(12) 

-0.38(11) 
0.125(65) 


Lat 

0.976(26) 
0.115(13) 
0.111(12) 
-0.52(14) 
0.211(90) 


TABLE L Energy of the two-pion system and the scattering phase shift 5. is the energy 
of the two free pions. Wn is reconstructed from the energy difference /SWn by Wn = AWn + W^. 
The invariant mass ^/s, the the momentum k and the phase shifts 5 calculated with the energy- 
momentum relations in the continuum (j24p are referred to as Cont, and those obtained with the 
relations on the lattice (j26p are referred to as Lat. The momentum is defined by k^ = s/A — m^. 
All values with the mass dimension are in units of the lattice spacing. 




vr 




P 




m 
E 


0.3567(24) 
0.6238(38) 




0.858(12) 
1.019(27) 




pIs/p^ 


Cont 

0.2620(45) 
0.956(16) 


Lat 

0.2799(50) 
1.021(18) 


Cont 
0.301(59) 
1.10(22) 


Lat 

0.358(74) 
1.31(27) 



TABLE IL Mass m and energy E of the pion and the p meson with the momentum p = pe^ = 
(27r/L)e3, extracted from the time correlation function. The momenta pcflf calculated from the 
dispersion relation in the continuum (I22p are referred to as Cont and that from the relation on 
the lattice (j23p as Lat. All values are in units of the lattice spacing. 
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